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Abstract. - We study the entanglement entropy of blocks of contiguous spins in non-periodic 
(quasi-periodic or more generally aperiodic) critical Heisenberg, XX and quantum Ising spin 
chains, e.g. in Fibonacci chains. For marginal and relevant aperiodic modulations, the entangle- 
ment entropy is found to be a logarithmic function of the block size with log-periodic oscillations. 
The effective central charge, Ceff, defined through the constant in front of the logarithm may de- 
pend on the ratio of couplings and can even exceed the corresponding value in the homogeneous 
system. In the strong modulation limit, the ground state is constructed by a renormalization 
group method and the limiting value of c B a is exactly calculated. Keeping the ratio of the block 
size and the system size constant, the entanglement entropy exhibits a scaling property, however, 
the corresponding scaling function may be nonanalytic. 



Introduction. — Recently, the entanglement proper- 
ties of strongly correlated systems have attracted great at- 
tention both in condensed matter physics [1] and in quan- 
tum information theory [2], for a review see [3]. Much 
work is devoted to homogeneous one-dimensional systems 
in which the entanglement entropy of L contiguous spins, 
defined as the von Neumann entropy of the density ma- 
trix of the block, Sl = — Trpi log 2 Pl, scales at the critical 
point as: 



Sl = - log 2 L 



(1) 



where c is the central charge of the associated conformal 
field theory and k is a non- universal constant. In the vicin- 
ity of the quantum critical point, L in Eq.(l) is replaced 
by the correlation length [1] . 

Inhomogeneities of different kinds (localized or extended 
defects, quenched disorder, etc.) are able to modify the 
local critical behavior of the system [4] and as a conse- 
quence, the scaling properties of Sl in Eq.(l) can also be 
changed. We expect that in case of weak perturbations 
caused by inhomogeneities the entropy scaling remains in- 
variant or is modified depending on the stability of the 
local critical behavior at the boundary of the block. For 
irrelevant perturbations the critical properties of the sys- 
tem at the boundary of the block remain unchanged and 



it is natural to assume that the same scaling law holds 
for the entanglement entropy as in the homogeneous sys- 
tem. On the contrary, for relevant perturbations the local 
critical behavior at the block boundary is governed by a 
new fixed point and consequently the scaling of the entan- 
glement entropy is expected to be modified. Finally, for 
marginal perturbations the local critical behavior is charac- 
terized by continuously varying local exponents, thus, also 
the prefactor of the logarithmic entropy scaling is presum- 
ably a continuously varying function of the strength of the 
inhomogeneity. 

This scenario has been checked [5] for a single defect 
coupling located at the boundary of the block in the spin- 
| XXZ chain defined by the Hamiltonian 



N 

Hxxz = ^ Ji{SfSf +l + SfSf 



(2) 



where the Sf's (a = x,y,z) are spin-i operators and 
Ji — J for i =/= L and Jl ^ J. Concerning the behav- 
ior of the entanglement entropy, this type of perturba- 
tion is found to be irrelevant in the ferromagnetic domain 
(A < 0) and relevant in the antiferromagnetic (AF) do- 
main (A > 0) in complete agreement with the local critical 
behavior of the system. More interestingly, the boundary 
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defect is a (truly) marginal perturbation in the XX chain 
(i.e. with A = 0), where the logarithmic scaling form in 
Eq. (1) is still valid, although c is replaced by a so called 
effective central charge, c e gf, which is found to depend con- 
tinuously on the strength of the defect [6]. We have ob- 
served [7] a similar marginal behavior with a continuously 
varying effective central charge in case of a single defect 
in the critical quantum Ising chain (QIC) defined by the 
Hamiltonian 

N N 

H I = -2j2JiS?Sf +1 -h'£S* (3) 

i=l i=l 

with Ji = J for i ^ L and Jl ^ J [8]. 

Quenched random disorder is a relevant perturbation 
both for the QIC and for the AF XX Z model [9]. An 
arbitrarily small random perturbation drives the system 
to an infinite randomness fixed point (IRFP), which can 
be studied in the framework of an asymptotically exact 
renormalization group (RG) method [10-12]. In accor- 
dance with this, it was found that the effective central 
charge jumps to the value characteristic for the IRFP for 
an arbitrarily weak random modulation [13-15]. 

The effect of non-periodic (quasi-periodic or more gen- 
erally aperiodic) disorder has many similarities to the ef- 
fect of quenched random perturbations. After the discov- 
ery of quasi-crystals [16], aperiodic systems have become 
the subject of intensive studies, both experimentally and 
theoretically [17,18]. Contrary to random disorder, these 
types of perturbations may also be irrelevant or marginal. 
The relevance of an aperiodic perturbation on the crit- 
ical behavior of a quantum chain is related to the sign 
of the cross-over exponent [19], <fi = 1 + v{w — 1), where 
v is the correlation length critical exponent of the pure 
chain and lo is the wandering exponent of the aperiodic 
model [20]. According to a heuristic criterion [19] the per- 
turbation is relevant if > 0, marginal if <j> = and 
irrelevant if (f> < 0. Later aperiodic QIC-s [21-23] and 
XX chains [24] were studied by an RG method and the 
low-energy properties were exactly calculated. For ape- 
riodic XXX chains, field-theoretical methods [25] and a 
variant of the strong disorder RG approach have been ap- 
plied [26-29] . For all models with marginal or relevant per- 
turbations anisotropic scaling behavior is observed which 
manifests itself in the scaling of the energy gap e ~ A~ z , 
where the dynamical exponent z is greater than one. 

In this paper, we study the entanglement properties of 
critical aperiodic quantum spin chains as a function of 
the strength of aperiodicity, measured by the ratio of the 
different couplings, r. We investigate systematically the 
effect of different type of non-periodic perturbations, such 
as irrelevant, relevant or marginal ones and examine the 
scaling form of the entropy in these critical, but not con- 
formally invariant systems. It is known that the spectrum 
of aperiodic systems have many unusual features [17] and 
we are interested in how these singularities are reflected 
in the entanglement properties. Spin chains both with 



continuous symmetry (XX and XXX models) and with 
discrete symmetry (QIC) are considered. 
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Table 1: Summary of sequences (Fib. - Fibonacci, Tripl. - 
tripling, P.D. - period doubling, Hier. - hierarchical) studied 
for different models (QIC, XX and XXX) in this paper. The 
relevance of the aperiodic perturbation at the fixed points (r = 
1 pure system, r — extreme aperiodic system) are indicated 
as I - irrelevant, M - marginal, R - relevant, whereas letter 
D refers to a noncritical dimerized ground state. The applied 
methods of investigations in the given range of r are also shown 
as: FF - free-fermion numerical calculation, RG - analytical 
renormalization group study. 

Methods. — In the case of QIC and XX models, which 
can be mapped to a system of free-fermions by means of 
standard techniques [30], we have performed large scale 
numerical calculations. Here, one determines first the re- 
stricted correlation matrix, Gf j, i, j = 1,2, . . . L, which is 
the corresponding L x L minor of the matrix defined in 
Eq. (2.32c) in Rcf. [30]. Then one determines the eigenval- 
ues vf (I — 1,2, . ..L) of the positive definite symmetric 
matrix [7] G L (G L ) T , where (G L ) T denotes the transpose 
of G L . In this representation, the entanglement entropy 
is given as a sum of binary entropies of the non-correlated 
fermionic modes [2,31]: 

S lW = - A ' lo S2( A i) " (1 " X i) lo g 2 (l - M), (4) 
l 

where A ; = (1 + ^)/ 2 - 

In the actual calculations, we considered blocks of size 
L which are given by finite approximants of the aperiodic 
sequence in order to get rid of log-periodic oscillations. 
The complete system then consists of N — 2L spins with 
periodic boundary conditions and the entropy is gener- 
ally averaged over the L different starting positions of the 
block. 

For the aperiodic sequences we used in this paper, we 
also consider the strong aperiodic modulation limit, r — > 
or 1/r — > 0, when the average entanglement entropy is 
studied by a variant of the strong disorder RG method. 
For the aperiodic AF XX and Heiscnbcrg chains, the 
ground state in this limit is of an aperiodic singlet form, 
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which is analogous to the random singlet phase of dis- 
ordered chains. The entanglement entropy is then given 
by the number of singlet bonds connecting the block with 
the rest of the system. For the strongly aperiodic QIC, the 
ground state is of an aperiodic embedded cluster form and 
the entropy is obtained by counting the clusters connecting 
the block with the rest of the system. The implementa- 
tion of the RG method is described in details where the 
specific problems are treated. The aperiodic sequences we 
study in this paper are selected in such a way to illustrate 
the general properties of the entropy in aperiodic quan- 
tum chains. The considered sequences are summarized in 
Table 1 , in which we have also given the stability of the 
two fixed points: the pure system's fixed point (r = 1) 
and the strongly aperiodic system's fixed point (r = 0). 
The applied methods used for the different models in the 
given range of the ratio r are also indicated. 

Fibonacci modulation. — Many basic features of en- 
tropy scaling can be seen for the Fibonacci modulation, 
which is irrelevant for the QIC (v = l,u> = — 1), marginal 
for the XX chain (y = 1, u = 0) and relevant for the XXX 
chain (y = 2/3, u> = 0). The Fibonacci sequence, which 
consists of two different letters a and b, is defined by the 
inflation rule: a — > ab and b — > a, so that we have by itera- 
tion: a, ab, aba, abaab, . . ., and the length of the sequence 
in the Z-th iteration is the Fibonacci number, F). In the 
quantum chains, the couplings take two values, J a = r and 
Jb — 1, depending on the underlying letter. To check the 
effect of an irrelevant perturbation for the QIC, we have 
calculated the size-dependence of the entropy for a strong 
modulation, r = 0.01, and plotted the results in the inset 
of the upper panel of Fig. 1 . We can see that the entropy 
has the same scaling form as in the homogeneous system 
with a central charge, c = 0.500(1). Hence this value is 
found indeed to be independent of the value of r. 

Repeating the same type of calculation for the Fibonacci 
XX model, the effective central charge is found to be r 
dependent, see the extrapolated values in the upper panel 
of Fig. 1 , which is in agreement with the marginal nature of 
the perturbation for this system. c c s(r) is monotonously 
decreasing with decreasing r < 1 and approaches a finite 
limiting value at r — with a correction of 0(r 2 ). 

We have calculated lim r ^n c e g(r) = c e ff(0) for the XX 
Fibonacci chain exactly by a variant of the strong disor- 
der RG approach [12]. For details of the application of 
the method, we refer to Refs. [26-29]. The basic idea is 
that two spins coupled by a very strong bond, Jb, form a 
maximally entangled singlet and the words aba and ababa 
are renormalized into the letters b' and a', representing 
the new effective couplings J' b ~ nr 2 and J' a w K 2 r 3 , re- 
spectively. For the XX model, the prefactor is k = 1. We 
note that this renormalization step corresponds to the re- 
versed triple application of the inflation transformation de- 
scribed earlier. The renormalized infinite chain is again a 
Fibonacci chain with J' a / ' J' b = r, so that one can repeat the 
transformation and finally the ground state of the system 
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Fig. 1: Upper panel: average effective central charge of the 
Fibonacci XX chain as a function of the coupling ratio, r. In 
the inset Sl vs. log 2 L is shown for the QIC at r = 0.01. Lower 
panel: the same for the tripling sequence. 

consists of a set of singlet pairs. The entanglement entropy 
in this case is given by the number of singlet bonds con- 
necting the block with the rest of the system [13]. The ef- 
fective central charge can be determined by calculating the 
difference AS = Sl — Sl> between the averages of entan- 
glement entropies belonging to blocks of length L = Fj +3 
and blocks of length U = F\. In order to do this, we no- 
tice that L/L' = p = t 3 for large I, where r = (1 + \/5)/2 
is the golden-mean ratio. The ratio of the length of the 
renormalized letters is X(a')/X(b') = t and the ratio of 
the density of the letters is given by /j,(a') / n(b') = t, thus 
the singlet bonds represented by letter b', cover a fraction 
n cov = 1/(t 2 + 1) of the chain. In this way, we obtain 

AS=^P-log 2 p = 2n cov , (5) 

and c cff (0) = 2/[(t 2 + 1) log 2 r] = 0.7962 which is in excel- 
lent agreement with the numerical findings (see the arrow 
in the upper panel of Fig.l). A more detailed derivation 
of Ccff based on the distribution of singlet lengths in the 
infinite system can be found in Ref. [32]. 

For the XXX-chain the Fibonacci modulation is a rele- 
vant perturbation and the fixed point at r = is strongly 
attractive. Indeed, the RG-procedure described for the 
XX-model leads to the same form of renormalized cou- 
plings but with a prefactor k = 1/2, thus the ratio 
J'a I 'J'b = r /2 < Ja/ Jb tends to zero even for a finite start- 
ing value of r. Consequently, the effective central charge is 
that given in Eq.(5) and it is independent of the coupling 
ratio r (provided 1). 

Tripling modulation. — The next sequence we con- 
sider is a tripling sequence, which is (for small pertur- 
bations) marginal both for the QIC and for the XX 
model (y = l,ui = 0) and relevant for the XXX model 
(u = 2/3, u) = 0). This sequence consists of three dif- 
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ferent letters, a, b and c and is denned by the inflation 
rule: a — > aba, b — > cbc and c — » abc. In the follow- 
ing, we consider Jb > J c > J a and use the parameter- 
ization: J c /Jh = J a /t/ c = r. For the QIC, the size- 
dependence of the entropy for r = 0.01 is shown in the 
inset of the lower panel of Fig.l, which is characterized 
by a central charge c e ff = 0.560(5). For the XX model, 
the shift in the effective central charge is even greater, 
see the extrapolated values in the lower panel of Fig.l. 
Surprisingly, this type of inhomogeneity enhances the en- 
tanglement as opposed to random and single defect per- 
turbations in these models. We mention that in some 
special chains of spins with many components also ran- 
dom disorder can enhance entanglement [33]. In the XX 
case, we have also calculated the small r limiting behav- 
ior by the RG method, in which singlets form over the 
Jb bonds and the new letters are defined by the deflation 
rules. The new couplings satisfy the marginality condi- 
tion: J' c / 1 J' b = J' a /J' c = r and in the aperiodic singlet 
ground state we have p = 3 and n cov — 1/3, such that we 
obtain from Eq.(5): c e ff(0) = 2/log 2 3 = 1.262 (see the 
arrow in the lower panel of Fig.l). For the XXX chain, 
the tripling perturbation plays a special role: a weak per- 
turbation at r = 1 is relevant but it is marginal at r — 0, 
where the same type of condition: J' b / J' c = J' c / ' J' a = r 
holds as for the XX chain. Consequently, the effective 
central charge of the tripling XXX chain is expected to be 
a continuously varying function of r at least for small r 
but discontinuous at r = 1. We note that the Fibonacci 
and the tripling XX (and XXX) chains belong to a class of 
aperiodic chains for which the ground state is an aperiodic 
singlet state in the limit r — > and the average entropy 
can be exactly calculated [32]. 

Period doubling modulation. In the rest of the 
paper, we focus on marginally non-periodic QIC-s for 
which even the first correction to the entanglement en- 
tropy at r = is analytically calculated, thus we can gain 
more insight into the varying nature of the effective cen- 
tral charge. First, we consider the period-doubling se- 
quence, which has the inflation rule a — > ab and b — > aa, 
such that we have by iteration: a, ab, abaa, abaaabab, 
abaaabababaaabaa, .... A small perturbation of this type 
is marginal for the QIC [y = l,w = 0), whereas in case 
of the XX and XXX chains, it drives the system into a 
non-critical phase with a dimerized ground state [24] . For 
the critical QIC, we use the parameterization: h c = 1, 
J a = r -1 / 3 and Jb = r 2 / 3 , when the dynamical exponent 
is given z = ln(r 2 / 3 + r~ 2 / 3 )/ln2 [21,22]. The numer- 
ically calculated effective central charge is shown in the 
inset of Fig. 2, which is found to be r-dependent. We have 
checked numerically that c c g(r) = c g(l/r). The limit- 
ing value is c e s(0) = 0, which is approached in a singular 
way as opposed to the smooth behavior of a single defect 
problem [6]. 

The limiting behavior for 1/r — > is studied in the 
frame of an RG approach [34], in which one transforma- 



tion step consists of two parts. In the first part, all strong 
Jb bonds are eliminated and two-site clusters with an ef- 
fective field, h w hi/ Jb = r~ 2 / 3 , are created. In the sec- 
ond part, all the strong transverse fields with h c = 1 are 
eliminated and new couplings J a w J 3 //! 2 = r _1 are cre- 
ated. Identifying Jb with J a = r -1 / 3 we have an aperi- 
odic sequence abbaaabbabb . . . which is in a period dou- 
bling form in terms of a' = a and b' = bb. During fur- 
ther renormalization steps the structure of this chain re- 
mains invariant. To calculate the entanglement entropy, 
we consider a chain with N — 4™ sites, in which there is 
a strong bond at L = N/2. By eliminating this bond, an 
effective cluster is formed, which connects the block with 
the rest of the system, such that in the limit 1/r — > 
Sl(N) = 1 and c e ff(0) = in accordance with the nu- 
merical results presented in Fig. 2. In the free-fermion 
representation, this result follows from Eq.(4) in which 
there is only one mixed non-correlated fermionic mode in 
the subsystem with a non-zero eigenvalue Ai = 1/2 (i.e. 
A; = 0, for I = 2,3, ... L). The perturbative correction 
to this term for 1/r <C 1 can be estimated as follows. 
For n = 2, i.e. for N — 16, we obtain the sequence 
abba after the first RG step, while by eliminating the two 
strong b bonds in the second RG step, a super-cluster is 
formed over the boundary of the block. The entanglement 
contribution due to this super-cluster is calculated per- 
turbatively leading to a new mixed fermionic mode with 
non-zero eigenvalue: A2 = r~ 2//3 /16. For n = 3, i.e. for 
N = 64, three RG steps can be performed, which results 
in a new super-super-cluster and a new mixed mode, so 
that there are then two non-zero sub-leading eigenvalues: 
A2 = A3 = r~ 2 / 3 /16. For a general n, we then have 
n embedded clusters and there are n — 1 non-zero sub- 
leading eigenvalues: A2 — A3 = . . . = A„ = r~ 2//3 /16. 
The entanglement entropy for small 1 /r indeed scales with 
n— 1 = log 2 iV/2 — 1 and the effective central charge is 
given by: c e ff(r) = r~ 2 / 3 (log 2 r/8 + const). The symmetry 
of the effective central charge c e s(r) = c e ff(l/r) observed 
numerically can also be seen in the perturbative treatment 
for small r. Hence, the function c e ff(r) is indeed singular 
at r = in accordance with the numerical results in the 
inset of Fig. 2. 

Hierarchical modulation. — The last sequence we 
consider in this paper is the hierarchical one in which the 
couplings are given by [35]: 

J t = Jr n , i = 2"(2m+l), n,m = 0,1,... (6) 

The hierarchical sequence is limit periodic [36] and can be 
generated through substitution with an infinite alphabet. 
For the QIC, the critical point is located at [8] h c = Jr 
and by setting J — 1/r, we have h c = 1. The excitation 
energy scales as e ~ N~ z , with a dynamical exponent 
z = ln(r+l/r)/ln2 [21,22]. (For the XX and XXX models 
this type of modulation drives the system to a dimerized 
phase.) 

In the numerical study we considered open Ising chains 
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of length N = 2™ and calculated the entanglement en- 
tropy of the half of the system, i.e. with L = N/2. As 
can be seen in Fig. 2, the entropy approaches a finite lim- 
iting value for any r / 1. Saturation of the entropy is 
found also at L = N/A and L — N/8 and expected to 
hold for any L — N/2 P , p — 1,2,..., since at these spe- 
cial positions there are very small couplings, which act 
in the infinite system as an effective cut. Repeating the 
calculation with a block of size L — [N/3]i n t, where [y}i n t 
denotes the integer part of y, we obtain the usual logarith- 
mic dependence as can be seen in Fig. 2. In this case c e s(r) 
is found to be a continuously varying function of r (see 
the extrapolated values in the inset of Fig. 2). We have 
thus the surprising conclusion that the scaling function 
S(x = L/N) = liniAr^oo Sl (N)/ log 2 N is non-analytic at 



x = 2-p, p = 1,2, 
mogeneous systems [1]. 



in contrast to the behavior in ho- 
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Fig. 2: Size dependence of the entanglement entropy of the 
hierarchical QIC for three different values r at L/N = 1/2 and 
for r = 0.2 at L/N = 1/3. Inset: extrapolated effective central 
charges as a function of r. Hierarchical chain at L/N = 1/3 
(0)> period-doubling chain at L/N = 1/2 (x). 

Next, we turn to a perturbative RG treatment for small 
r and focus on the block composed of the first L spins 
at the left boundary but now we consider the thermody- 
namic limit N — > oo. As for the period doubling modula- 
tion, we eliminate the strongest bonds, which are at odd 
positions (i.e. with n = 0), and create renormalized two- 
spin clusters subjected to an effective transverse field of 
h w h 2 / J = r. The renormalized chain has now N/2 sites 
with the same hierarchical structure of couplings. In the 
repeated use of the transformation, bonds with n = 1 
(then with n — 2, etc) are decimated out and super- 
clusters with a hierarchical structure are created. Then for 
any value of L, the block is connected to the environment 
through one super-cluster, thus the entanglement entropy 
is Sl = 1 and c e ff(0) — in accordance with the results 
presented in the inset of Fig. 2. This result can also be de- 
rived using the free-fermion representation of the system, 
when in the limit r — > 0, the density matrices of the non- 



correlated fermionic modes have eigenvalues (Aj, 1 — Aj) 
given by Ai = 1/2 and A; = 0, for I = 2,3,... L, thus 
S L = 1. 

For a small but finite r, there are perturbative correc- 
tions, which depend on the value of L. In the simplest case 
L = 3, these corrections are due to correlations between 
the spin-clusters (1,2) and (3,4). In the non-correlated 
free- fermionic description, the eigenvalues are up to 0(r 2 ): 
Ai = l/2 + 0(r 2 ), A 2 = r 2 /4 and A; = 0, for I = 3, . . . L. 
For L = 5 there are two perturbative cluster-cluster con- 
tributions which will result in A 2 = 2r 2 /4, whereas the 
other eigenvalues remain unchanged, at least up to 0(r 2 ). 
Generally, for L K = 2 K + 1, k = 0, 1, . . . there are k per- 
turbative contributions and in the limit nr 2 <C 1 we have 
A 2 = nr 2 /A. (cf. the form of perturbative cluster correc- 
tion in the case of the period doubling sequence.) Thus 
we obtain for the entanglement entropy: 



?2~ + l 
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O(kV) 



(7) 

For a general value of L = 2 n (2m + 1), the number of 
independent perturbative (cluster) contributions, k, does 
not depend on n and given by k — [log 2 ra]j„ t + 1. In 
this case the entanglement entropy in the leading order is 
given in Eq.(7). We thus conclude, that Sl{t) has log- 
periodic oscillations in L and the leading L dependence 
is faster than logarithmic in the range (log 2 L)r 2 <C 1, 
although for (log 2 L)r 2 3> 1, a logarithmic L dependence 
is expected. 

Summary. — Summing up, we have studied the en- 
tanglement properties of critical quantum spin chains, 
such as the antiferromagnetic Heisenberg and XX mod- 
els, as well as the quantum Ising chain in the presence of 
aperiodic modulations of the couplings. As far as the crit- 
ical properties of the systems are concerned, an aperiodic 
inhomogeneity can be irrelevant, marginal or relevant, de- 
pending on the model and on the fluctuation properties 
of the aperiodic sequence. An irrelevant aperiodicity al- 
ters neither the logarithmic scaling of the entropy nor the 
constant in front of the logarithm, i.e. the growth of the 
entropy is determined by the central charge c of the pure 
system. For marginal and relevant modulations, the av- 
erage entropy is still scaling logarithmically with the size 
of the block, however, the prefactor is modified compared 
to that of the homogeneous system i.e. c e ff ^ c and in 
addition to this, log-periodic oscillations appear. Inter- 
estingly, c ff can exceed the corresponding value in the 
homogeneous system and for marginal perturbations, it 
is a continuous function of the coupling ratio r. In the 
r — > limit, where the ground state of several aperiodic 
sequences can be exactly constructed, the average effective 
central charge was analytically calculated. 
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